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Abstract 

Keldysh nonequilibrium Green's function method is utilized to study theoretically the spin polar- 
ized transport through a graphene spin valve irradiated by a monochromatic laser field. It is found 
that the bias dependence of the differential conductance exhibits successive peaks corresponding to 
the resonant tunneling through the photon- assisted sidebands. The multi photon processes origi- 
nate from the combined effects of the radiation field and the graphene tunneling properties, and 
are shown to be substantially suppressed in a graphene spin valve which results in a decrease of the 
differential conductance for a high bias voltage. We also discussed the appearance of a dynamical 
gap around zero bias due to the radiation field. The gap width can be tuned by changing the 
radiation electric field strength and the frequency. This leads to a shift of the resonant peaks in 
the differential conductance. We also demonstrate numerically the dependencies of the radiation 
and spin valve effects on the parameters of the external fields and those of the electrodes. We 
find that the combined effects of the radiation field, the graphene, and the spin valve properties 
bring about an oscillatory behavior in the tunnel magnetoresistance(TMR), and this oscillatory 
amplitude can be changed by scanning the radiation field strength and/or the frequency. 

PACS numbers: 85.75.-d, 81.05.Uw, 75.47.-m 
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I. INTRODUCTION 



Laser fields have been demonstrated to be a powerful tool for the exploration and modifi- 
cation of the materials properties. For example, in conventional semiconducting nanostruc- 
tures a variety of applications have been established such as the various forms of the photo 
electronic devices (e.g., the radiation-controlled field-effect transistors, photodiodes, and the 
light-emitting diodes 1-5 ). On the other hand, graphene with its many unusual physical prop- 
erties (to mention but a few, the half integer quantum Hall effect^, the Klein tunneling^ and 
the conductance properties 8 ) offers a new platform to explore radiation effects. The influ- 
ence of electromagnetic fields on graphene has been recently a subject of intense research. 
Particularly interesting examples are the linear response to and the frequency dependence 
of the conductivity^ - — , the photon-assisted transportr^ 1 ^, the microwave and far-infrared 
response^ - —, the plasmon spectrum^ - —, as well as the nonlinear response to the electro- 
magnetic radiation 20-23 . The potential applications of graphene in terahertz electronics were 
pointed out in Refs.—. 

The resonant interaction between graphene and an electromagnetic field may open a 
dynamical gap in the quasiparticle spectrum of graphene 2 ^, which leads to a strong sup- 
pression of the quasiparticle transmission through a graphene p-n junction. However, a 
directed current without applying any dc bias voltage can be generated in certain conditions 
as a result of inelastic quasiparticle tunneling assisted by one- or two-photon absorptions 2 ^. 
Another way to generate a direct current is to use shaped, timely tuned electromagnetic 
pulses^ which also may generate a valley current. Therefore, it seems possible to control 
the transport properties of diverse graphene tunneling structures by the variation of the 
strength and the frequency of the external radiation fields. On the other hand, the potential 
of graphene for spin- dependent transport (spintronic) applications is well documented by 
now2£~— . Motivated by these facts, in this work, we study theoretically the spin polarized 
transport through a graphene-spin-valve device in the presence of a monochromatic laser 
field. The method is based on the standard Keldysh nonequilibrium Green's function ap- 
proach, as described in Refs.— and^I. We find that additional peaks emerge in the bias 
dependence of the differential conductance that reflect the resonant tunneling through the 
photon-induced sidebands. The strong suppression of the multiple photon processes stem- 
ming from the combined effect of the radiation field and the tunneling through the graphene 
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leads to a strong decay of the differential conductance at a high bias voltage. The resonant 
interaction of the quasiparticle in graphene with the radiation field results in a dynamical 
gap in the quasiparticle spectrum. When the bias voltage lies inside this gap, the differential 
conductance displays a zero value region situated symmetrically around the zero bias. The 
width of this region can be tuned by changing the radiation strength and the frequency, 
which also cause a shift of the resonant peaks in the differential conductance. The TMR 
versus the bias voltage exhibits an oscillatory behavior, and its oscillatory amplitude can be 
controlled by the radiation field strength and/or the frequency. 



II. THEORETICAL MODEL 



We consider a spin valve device consisting of an extended graphene sheet (that defines 
the x — y plane) contacted by two ferromagnetic electrodes. A gate voltage V g applied 
on graphene shifts the Dirac point away from the zero energy. The two electrodes are 
voltage-biased with respect to each other with a bias V. The electrical current flows in 
the x direction. Additionally, we assume that a laser field is homogenously irradiating the 
structure. Within the metallic electrodes we assume the field to be shielded and ignore thus 
its effect on the electrodes. The laser field is monochromatic and is linearly polarized along 
the y direction. The low energy graphene Hamiltonian around the Dirac points has the 
minimal coupling form 

Hq = vfct • k — evFA(t)(7y, (1) 

where k is the operator of the electron momentum in the graphene plane, e is the electron 
charge, cr is the vector build out of the Pauli matrices in the sublattice space, and A is 
the laser's vector potential. In the tight binding description^, vp = ^at g with t g being 
the nearest neighbor hopping energy, and a being the carbon-carbon distance. The vector 
potential A(t) is taken as 

E 

A(t) = — sinw t, (2) 

where Uq is the frequency of the radiation field, and E is its amplitude. Diagonalizing the 
Hamiltonian (TjQ) in the absence of the radiation field, one finds the eigenvalues 

e sk = sv F \k\ +V g (3) 
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with s = ± denoting the band index, and the eigenstates 

1 



iM') = ^f* k X, (4) 



where Q is the volume of the system, and 

< = ^(j<*>)' tm0(k) = l' (5) 

We introduce the field operators 

# T (r, t) = J2 ^skr(t)ip s k(r), (6) 

where a skr is the usual annihilation operator for an electron in the band s, with the momen- 
tum k and the spin r, and then express the Hamiltonian ([T]) in the second quantized form 

as 

H G = Yl e s k<^ kr a sk T - ev F A(t) Y dss'^lkr^'kr, (7) 

ksr kss'r 

where 

d ss/ ^= l -(se-^-s'e^). (8) 

Applying the rotating wave approximation 39 , we neglect the energy nonconserving terms. 
Accounting for the coupling between the graphene and the two ferromagnetic electrodes we 
write for the Hamiltonian of the complete tunnel junction in the presence of the laser 

H = H G + H L + H R + H T , (9) 



where 



Hq — Y e sk«lkr a skr 
ksr 

-ev F £ <U*e-^+_, k < kr a_ kT (10) 

kr L 

+A e^°*rf_ +ik a t _ kT a +kT } 



with A = —i-^-, and 



H\ = y^gqATCq Ar C q Ar, X = L, R, (11) 
qr 

H T = (^qAkC qAr a sk r + #.c.) . (12) 
qkArs 

Eqs. (TTTT) and (TT2T) describe respectively the A electrode and the coupling between the 

graphene and the electrodes. e n \ T is the single electron energy, c qAr (c q A T ) is the usual 



creation (annihilation) operator for an electron with the momentum q and the spin r in the 
A electrode. 

By introducing a unitary transformation 



U = exp 



,u t 
—% 



■) y^( Q +kr Q +kr - aL kr a-kr) 
kr 



(13) 



we redefine the Hamiltonian of the system in the rotating reference as 



H = U^HU + i^ r U 



— ^2 eskQskT^kT + Yl A (a+ kr Ct-kr + «^kr a +kr 
kr kr V 

+ ££qAr< Ar C qAr + £ fekW<Ar^kr + ^.C 



(14) 



where 



and 



qAr qkAsr 



Csk — Csk — SUo/2. 

In the calculation of Eq. (ll4jl . we have assumed (as in Refs.— >2S»ai) that the most impor- 
tant contributions stem from an almost one dimensional electron motion (k x 3> k y ) in the 
interaction with the radiation field. 

The electric current of the system can be calculated by the time evolution of the occu- 
pation number operator of the left electrode, 

/ = e(Xf L ) = l ^([H,M L }), where M L = ^ < Lr c qLr . (15) 

qr 

Using the nonequilibrium Green's function method, Eq. f|T5|) can be further expressed as 
I = -{E/*JSM[(^(t.l4)-5^(t l ti))t(e) 

16 

+^(t,ti)]rL}e-^-*)e~^ L e^, 



where Tr is the trace in the spin space, f\(e) is the Fermi distribution, Qggi(t,ti) = 
X^kskvCMi)' an d Gfgr(t,ti) = X^k.skv(Mi) are 2 x 2 matrices denoting the retarded 

kk' ' kk' 

(advanced) Green's function and the lesser Green's function, respectively. In the calculation 
of Eq. ffTBT) . we assume that the dominant contributions to the tunneling stem from the 
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electrons near the Fermi level, and hence the linewidth function to be independent of k. 
Thus, we have 



rl o 
o ri 



(17) 



with = 27r Tq Ak T q Ak'5(£ — £ q Ar)- in order to solve Eq. (|T6|) . we need to calculate the 
q 

Green's functions Q T J s i ,r (t,t') and Q T J s i ' (t,t r ). Using the equation of motion method, we get 

+ J dt x dt 2 g^(t - t 1 )E^(t 1) t 2 )^'^(t 2 ,t'), 



where 



d,B , j-i-l .-stunt .s'wnt' 



K^t') = J ^e-^e^e^ (19) 

with Sq r = — f (F£ + rjj), and g TT '' r (t) is the retarded Green's function of graphene without 
the coupling of the electrodes, and can be obtained by a straightforward calculation. The 
detailed expressions are given in the Appendix. Q TT '<(*,£') is related to Q TT ' r {t,t') through 
the Keldysh equation 

g TT '<(t,t>) = fdhdtzG^fatjZ^i^tJG^fat'), (20) 

where 



K' s >(t,t') = iI^[rih{e) + Vlf R {e)\ 

■ (J. y7\ ■ SUJQt .S UQt V / 

xe~ K1 >e l ^~e~ % 2 . 



Notice that the Green's functions in Eqs. (|T8j) and (|20|) do not depend only on the difference 
of the two time variables, thus one should take a generalized Fourier expansion as^ 

g(t, = ^2 / dee- iet e i{£+nuJo/2)t 'g(e } e + nw /2). (22) 

n 

Hereafter we shall use the simple notation 

g n m{£) =g{e + nuo/2, e + mu /2). 
Evidently, the different Fourier components satisfy the relation 

g n m{e) = g n -mfl{e + muo/2). 
From Eq. f fT8|) . the Fourier components of the Green's function can be expressed as 



=?rn,n "mn 



(23) 



v. 



m,m—2 



m,m+2 



where e m = g%£E ' , and 

i/?Ttm;H — IW^O u 

9mm; ( £ )^0 

fi'mm;++( £: )^V 
fi'mm; — \-( £ )^0 

with g T ^nss'( £ ) = 3mngll> r + m Uo/2). Note that Eq. (|23|) is formally equivalent to the ones 
describing the motion of electrons in a tight-binding linear chain with site energies and the 
nearest-neighbor coupling. Its solution can be derived by using the conventional recursive 
technique^. 

Carrying out the generalized Fourier transformation to Eq. fl20|) . and substituting it in 
Eq. (|T6l) together with Eq. (|23|) . we finally obtain the time-averaged current (under time 
average, the displacement current for the central region is zero since it is given by the charge 
accumulation in the central region over the time period and tends to zero for large time 
durational) 

(i) = -fE/£{^oT'( £ )-^ T ;n a (^ s K,o 

TTl 

+^ < ^)s T L ;,o(^)} ) 



(24) 



where 



Gl^{e) = £ £ ^;(^)S- n < n2 (e 2 )^ 2 ;;; a (e 2 ), (25) 

T\OL n\Tl2 



S mn fx(s+^^) 5 n , m+2 f x (E + ^^)\ (26) 
S n>m -2fx(e + Kmfx{e + ^^) 




^imn( e ) = lT X I . • (27) 



In Eq. (1241) . we further set the symmetrical voltage division: fx L . R = E F ± l/2eV, and put 
Ep = in the numerical calculations. The expression (124")) is the central result of this paper, 
and allows one to describe the spin-polarized transport in a nonlinear response regime, i.e., 
high bias voltage or high strength of the radiation fields. Please, note the difference to our 
previous study^ in which we studied the influence of a time-dependent chemical potential 
(realized as an oscillating bias voltage). In the present study the radiation field, coupled 
to the carrier as a vector potential, induces a dynamical gap in the quasiparticle spectrum. 



Thus, in turn leads to a marked changed in the transport properties and TMR, as shown 
below. The TMR is deduced according to the conventional definition 



where l(0(n)) is the time-averaged current flowing through the system in the parallel (an- 
tiparallel) configuration. 

III. NUMERICAL ANALYSIS 

Based on the above analytical expressions, we present the numerical calculations for the 
spin polarized transport through the system. Firstly, we assume the linewidth function 
r^(e) to be independent of the energy within the wide band approximation, and the two 
electrodes are made of the same material. Introducing the degree of the spin polarizations 
of the left and the right electrodes Pl = Pr = P, we can write = rj^ = r (l ± p) where 
r describes the coupling between the graphene and the electrodes in the absence of an 
internal magnetization. We set the temperature to zero and t g ~ 2.8 eV in our calculation. 
In addition, we choose the frequency huo < 0.06t 9 , and the voltage eV < lAt g so that the 
Dirac approximation remains valid for such energies^. 

FigJT](a) shows the bias dependence of the differential conductance G = d{I)/dV for 
different radiation strengths in the parallel configuration of the two electrodes. In the ab- 
sence of the radiation field, the differential conductance exhibits a strong dip structure (see 
the Fig{T](a)). The graphene spin valve device is sensitive to the electromagnetic radiation. 
When the radiation field is applied, some characteristic features are presented explicitly in 
the transport. Firstly, close to the zero bias voltage, a salient zero value region appears 
symmetrically in the differential conductance, and the width of the region becomes large 
with increasing the radiation strength. This is because the resonant interaction between 
the graphene and the external radiation field opens a dynamical gap in the spectrum of the 
quasiparticle in the graphene. When the bias voltage lies inside this gap, the graphene can 
be viewed as an insulator, thus the electron tunneling through the graphene is suppressed. 
Additionally, the combined effects of the radiation field and the graphene tunneling prop- 
erties results in multi photon processes (see Eq. ( |23|) ). which are manifested a successive 
peaks in the differential conductance as a function of the bias voltage. The peaks close to 



TMR 




(28) 



1(0) 
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and far away from the zero bias point are suppressed. This is due to the nature of the 
insulator-like properties in the dynamical gap (for peaks close to eV = 0), and (for those far 
away from eV = 0) due to the lower amplitude of multi photon processes with an increasing 
number of exchanged photons at a fixed field strength. In general, with the increase of 
the radiation field strength, the resonant peaks are weakened (see FigJT^b)). This behavior 
is different from the case of when an oscillating gate voltage is applied, in which case the 
resonant peaks increase uniformly with increasing the ac gate voltage amplitude^. The 
reason for this difference is that the sidebands for an oscillating gate voltage are produced 
due to the modulation of each of the quasiparticle levels by the gate voltage, while those 
for the present case are caused by the electron scattering between the bands. When the 
electrons are injected from the electrode to the graphene, the radiation field results in mul- 
tiple electron scattering events. Each scattering process provides a sideband. When the bias 
windows cross this sideband, the resonant tunneling occurs leading to the appearance of 
the resonant peak. When the radiation strength is enhanced, the dynamical gap increases. 
Thus the sidebands fall into the gap at some stage, and the scattering probability of the 
electrons is suppressed due to the linear density of state (DOS) in graphene which explains 
the decrease of the peaks. Another pronounced feature is that the positions of the peaks 
shift in the direction away from the origin with increasing the radiation strength, which also 
differs from the case of a gate voltage 13 , where the positions of the peaks are fixed. This 
is due to the increase of the dynamical gap lifting the photon-assisted sidebands. Fig. [flc) 
shows the bias dependence of the differential conductance for different frequencies uq in the 
parallel configuration of the electrodes. With increasing the frequency, the dynamical gap 
diminishes leading to a decrease of the width of the zero value region and to a shift of the 
peaks towards the origin. Especially, one can find that the amplitude of the peaks become 
larger when the frequency grows. The reason is that the frequency increase lifts the photon- 
induced sidebands, which causes the increase of the DOS for the sidebands, thus enhancing 
the electrons tunneling through the graphene sheet. 

Figf2](a) shows the bias dependence of the differential conductance for the different polar- 
izations p in the parallel electrodes configuration. In the differential conductance the peaks 
are viewed as stemming from the resonant contributions of the sidebands. And the other part 
may be attributed to the off-resonant tunneling. It is observed that the latter contribution 
is almost independent of p; while the former exhibits a sizable dependence on the polariza- 
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tion. This can be understood roughly in that the sidebands become more spin- selective with 
increasing the polarization, thus gradually turning off the spin-flip channel (mixing channel) 
on the sidebands. This decreases the differential conductance for the resonant states. For 
very weak peaks, the contribution from the background off-resonant tunneling may com- 
pete with those from the resonant states, making the total change vague. Figf2^b) shows 
the bias dependence of the differential conductance for different gate voltages. It is found 
that the oscillation peaks and the zero value region shift towards the positive direction of 
the bias voltage since the positions of the Dirac points vary due to the gate voltage. The 
additional interesting phenomenon is that the amplitude of the oscillation peaks increases 
with increasing the gate voltage. This is because the gate voltage lifts the quasiparticle 
level of the central graphene in the bias windows, the linear DOS of graphene enhances the 
tunneling probability from the ferromagnetic electrode to graphene. 

The differential conductance as a function of the radiation strength for different radiation 
frequencies in the parallel configuration of the electrodes' magnetizations is shown in Fig. 
Ela). The differential conductance with varying the radiation field strength exhibits several 
resonant peaks. This can be traced back to the fact that, with a fixed bias window and 
a fixed frequency, the dynamical gap in the quasiparticle spectrum of graphene becomes 
larger for stronger radiation field strength, which pushes the sidebands upwards. When the 
sidebands cross the bias windows, the resonant transport occurs, which is then manifested 
in the conductance peaks. When the radiation field strength is large enough, the differential 
conductance tends to a constant value which almost does not vary with the frequency, hinting 
on the conclusion that graphene can be identified as a pure resistance in this situation. 

Fig. Elb) shows the dependence of the differential conductance on the radiation frequency 
cuo for the different radiation strength in the parallel configuration of the electrodes. It is 
clearly seen that when uq approaches zero (note the magnitude of the field, i.e. E Q , is finite), 
the differential conductance diminishes owing to the large dynamical gap. With increasing 
the frequency, the dynamical gap shrinks, which causes lowering in energy of the photon- 
induced sidebands. Therefore, this inverse process (to Fig. EI a)) causes the resonant peaks 
as well in the spectrum in the bias windows. It is found further that the interval between the 
peaks is not uniform since the radiation frequency changes the dynamical gap in a nonlinear 
way. Moreover, increasing the radiation field strength enlarges the dynamical gap, and 
moves the resonant peaks in the differential conductance away from the origin. 
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The bias dependence of the TMR for the different radiation strengths and frequencies is 
shown in Fig. HI The TMR as a function of the bias voltage exhibits an oscillatory behavior. 
We assign this behavior to the combined effects of the radiation field, the graphene and spin 
valve properties. In order to clarify this point more clearly, We plot the bias dependence 
of TMR for the different polarization in FigJSJ One can observe that the TMR changes 
in a non-linear manner: the TMR values in the off-resonant regions become larger than 
those at resonant points. This is because graphene remains the insulator-like in the off- 
resonant tunneling, the ballistic spin tunneling enhances the TMR. In addition, at the 
high bias voltage, the oscillation peaks decay sharply due to the suppression of the multi 
photon processes in graphene. With increasing the radiation strength (or frequency), the 
oscillation is enhanced (or weakened) due to the interplay between the spin valve effect and 
the dynamical gap (see Fig. 0J. 



IV. SUMMARY 



In conclusion, we studied theoretically the spin polarized transport through a graphene 
spin valve device assisted by a linearly polarized, monochromatic laser field and in the 
presence of a dc bias and a gate voltage. The method is based on the standard Keldysh 
nonequilibrium Green's function approach. We find that the bias dependence of the dif- 
ferential conductance exhibits additional peaks due to the resonant transport through the 
photon-assisted sidebands. The resonant interaction of the quasiparticle in the graphene 
sheet with the radiation field causes a dynamical gap in the quasiparticle spectrum. When 
the bias voltage lies inside this gap region, the differential conductance displays a zero value 
region situated symmetrically around the zero bias. The value of the dynamical gap de- 
pends on the strength and the frequency of the external radiation field, thus the width of 
this zero value region in the differential conductance can be tuned by changing the radiation 
strength and the frequency, which also causes the shift of the resonant peaks in the differ- 
ential conductance. We explored the behavior of the peaks in the differential conductance 
with varying the radiation field strength and its frequency and exposed the dependence on 
the spin polarization of the ferromagnetic electrodes. We also demonstrated that the com- 
bined effects of the radiation field, the graphene and the spin valve properties result in an 
oscillatory behavior in the TMR. This oscillatory amplitude can be varied by changing the 



11 



parameters of the radiation field. 
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Appendix A 

In this appendix, we present the analytical results of the Green's function g T s l,' r {t) for 
the graphene without the coupling to the electrodes. From the Hamiltonian H G , we find by 
using the equation of motion method that 



9 T s 7(t) = J 



(Al) 



where g TT ' r (e) = Y / 9k' r ( e ) with 

k 



£-e_ k 



TT,r ( \ I ( £ - e+k )(e-(E_ k )-A 2 (e-e +k )(e-e_ k )-A 2 



(A2) 



(e-? +k )(e-?_ k )-A 2 (e^ +k )(e^_ k )-A 2 

To solve glli' r (s), we need to convert the summation over k into a integral in two dimensional 
momentum space. After a straightforward calculation, we have for \e\ < A, 



TT,r 

9u 



irv'i, { 



e-V„ 



In 



{D-V g -uj /2) 2 +A 2 -{e-V g ) 2 



c 2 /4+A 2 -( £ -V 9 ) 2 



- v/A 2 - (e-ig 2 [arctan 

— arctan 



1/2 



D-V g -L> a /2 

^A 2 ~{e-V g ) 2 

+ D-V g }, 



TT,r 
922 



1 r (e-V g ) 



2ttv' 



In 



(D-V g -wo/2) 2 +A 2 -(e-V g ) 2 



uj 2 /4+A 2 -(e-V 9 ) 2 



+ v /A 2 -(e- V^[arctan ^g=| 

D + V g }, 



arctan 



-^0/2 I 



TT,r TT,r 
9\2 =921 



4nv 



Kin 



(D-V g -u] Q /2) 2 +A 2 -(e-V g ) 2 
u 2 /4+A 2 -( £ -V g ) 2 



(A3) 



(A4) 



(A5) 
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For \e\ > A, 



tt ,r 

9u 



4lTV 



(D-V, 



/2 )2-( e -y g )2+A 2 



V(e - V g ) 2 - A 2 

(D-\/ g -^/2 +v /( £ -V 9 )2-A2)(o, /2 +v /( £ -y 9 )2-A2) 



(D-V 9 -cu /2- v /(£-V 9 ) 2 -A2)( W0 /2- v /( £ -y 9 )2-A2) 



x In 



-2V g ] + ^sgn(e - V g ){£(A < \e - V g | < + A 2 ) 



+ 2L> 



x(e-V s -v/( £ -^) 2 -A 2 ) 

-0(A < |e - < v /( J D-K ; -f) 2 + A 2 ) 

*(e-V g +y/(e-V g y-A*)}, 



tt ,r 
922 



A-KV 



Szfe-V g )hi 



(D-V g -u Q /2) 2 -(e-V g y 2 +A 2 
( £ -V 9 )2-A2- w 2/4 



x In 



-2D 



+ y / (e-V 9 ) 2 -^ 2 

(g-y 9 -o;o/2+ v ^y 9 )2-A2)(a; /2+ v /( £ -\/ 9 )2-A2) 

(D-y 9 -o,o/2- v /( £ -y 9 )2-A2)( t u () /2- v /( £ -y 9 )2-A2) 
+2F 9 ] + ^sgn(e - V g ){£(A < |( e - V g )| < y/* + A 2 ) 
x(e-V g + y/(e-V g )*-A*) 



-9{A<\e-V g \< ^{D-V g - 
x(e-V g -^(e-V g y-A*)} : 



<^0 \2 
2 > 



A 2 ) 



9ll ' T = 9ll ' r 



(A6) 



(A7) 



47 ™f- v /( £ -V 9 ) 2 -A 2 



(g-V 9 -^o/2+ v /( £ -\/ g )2-A2)(u;o/2+ v /( £ -V 9 )2-A2) 

(D-V g -w /2-y/( e -V 3 )*-A*)(u, /2-y/( e -V g )*-A*) 



X In 



+i Vsgn(5 - V g ){£[A < |e - V g | < V? + A 2 ] 



V(£ ,^)2, A 2 - ^ < l g - V '\ < -V 9 - f ) 2 + A 2 ] 



V /(^V 9 )2-A2- 

D is a high-energy cutoff of the graphene bandwidth. 



(A8) 
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FIG. 1: (color online) The bias dependence of the differential conductance G for different radiation 
strength Eq at fluo = 0.04t g (a), and for different radiation frequency ojq at Eq = 1200V/ cm (b). 
The other parameters are taken as p = 0.4, To = 0.05t g , V g = and D = 3t g . The inset in (a) 
shows the differential conductance in the absence of the radiation field. 
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FIG. 2: (color online)The differential conductance G versus the bias voltage (a) for different po- 
larizations p at wo = 0.03t g and Eq = 1200V/ cm,V g = 0, and (b) for different gate voltages V g at 
ujq = 0.04t ff and Eq = 1200 V/ cm,p = 0.4 . The other parameters are taken the same as those of 
Figffl 
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FIG. 3: (color online)The differential conductance G as a function of the radiation field strength 
for different frequencies (a), and as a function of the frequency for different radiation strengths (b). 
In both cases we consider the parallel configuration of the electrodes magnetizations at eV = 0.5t g . 
The other parameters are taken the same as those of Fig{TJ 
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FIG. 5: (color online) The TMR versus the bias voltage for a different polarization p at Hloq = 0.06i 9 
and Eq = 400V/ cm. The other parameters are taken the same as those of FigfTJ 
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